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1 Introduction 

According to the compilation of the Particle Data Group 

(PDG) 2002 [1], the absohitc vahics of the entries in the 
first row of the Cabibbo-Kobayashi-Maskawa (CKM) mix- 
ing matrix are given by 

\Vud\= 0.9734 ±0.0008, 
\Vus\ = 0.2196 ±0.0026, 
IKbl = 0.0036 ±0.0010, (1.1) 

which implies a 2.2 a deviation from unitarity: 

\Vud? + IKsl' + IKbl' - 1 = -0.0042 ± 0.0019. (1.2) 

The value for \ Vud\ in (1-1) has been extracted from super- 
allowed Fermi transitions of several 0^ nuclei and neutron 
beta decay, whereas the number for \Vus \ is based on more 
than thirty- year-old Kg^ data. 

The situation has changed dramatically with the out- 
come of a new high statistics mcasiirement of the A'i^ 
branching ratio by the E865 Collaboration at Brookhaven 
[2]. Their analysis of more than 70,000 events yielded 
a branching ratio which was about 2.3(7 larger than the 
current PDG value. As a consequence, the value of \Vus\ 
based on the new experimental result does not indicate 
any significant deviation from unitarity. Moreover, besides 
indicating a sharp disagreement between new and old K^^ 
data, the new result implies an inconsistency between if 
and ifgg data. 

The current experimental information on the decay 
mode of the neutral kaon is indeed very unsatisfactory. 
The two numbers given by PDG 2002 [1], 



(7.50 ±0.08) X 10^ s- 
= (7.7 ±0.5) X 10S-\ 



(1.3) 



differ considerably depending on the procedure for the 
treatment of data. The first value in (1.3) was obtained 
from a constrained fit using all significant measured 

branching ratios, the second one is a weighted average of 
measurements of the /sTgg ratio only. Apparently, the rate 
obtained from the fit is completely driven by input differ- 
ent from the actual measurements. In particular the error 
on the "fitted" value does not reflect at all the experimen- 
tal accuracy (the experiments were made in the sixties and 
early seventies) but rather the constraints from the global 
fit. 

Presently, new independent ffes decay measurements 
are in progress (CMD2, NA48, KLOE) and should help to 
clarify the experimental situation. 

In this paper, we present a detailed numerical analysis 
of the radiative corrections to the ii'^g Dalitz plot distribu- 
tion. We discuss possible strategies to extract \Vus\ from 
the experimental data and we propose a rather powerful 
consistency check of K'^^ and K'^^ measurements. 

This work is based on our previous calculation [3] of 
the Ki'i decays to 0{p^ . (md — mu)p'^ , e^p^) in chiral per- 
turbation theory with virtual photons and leptons [4] . Af- 
ter a brief review of the main kinematic features of i^es 



decays and the structure of radiative corrections (Sect. 2), 
we recall the structure of the form factors relevant for i^eS 
decays including a discussion of the recent results on the 
contributions of order in the chiral expansion in Sect. 3. 
Real photon emission in the K^^ case is discussed in Sect. 
4. In Sect. 5 we illustrate our general considerations by a 
numerical study of the K^^ decay and the description of a 
procedure to extract the CKM matrix element \Vus\ from 
experimental data. The impact of the E865 experiment 
on the determination of \Vus\ from iiT^ data is discussed 
in Sect. 6. A specific strategy for a combined analysis of 
ifgg and data is proposed in Sect. 7. Our conclusions 
are summarized in Sect. 8, and three Appendices contain 
some technical material related to the calculation of loop 
contributions and real photon radiation. 



2 Kinematics and radiative corrections 

The generic ii'es decay 

K{pk) ^ 7r(p^) e+(pe) MP'^) (2-1) 

can be described by a single form factor (usually denoted 
by /+). A second form factor^, being also present in princi- 
ple, enters only together with the tiny quantity mg/Af|- ~ 
10~^ in the formula for the Dalitz plot density. Therefore, 
these contributions are utterly negligible and the invariant 
amplitude (in the absence of radiative corrections) can be 
simplified to 

M = ^ V:, C ]f it) {pK + P.)„ (2.2) 



where 

= u{p,) (1 - 75) v{pe) 
denotes the weak leptonic current, and 

1 for 



C = 



I/V2 for if+3 • 



(2.3) 



(2.4) 



The form factor depends on the single kinematical variable 
t = [pk — P-kY ^'^d the superscript (0) indicates the limit 
e = 0. 

The spin-averaged decay distribution p{y,z) for K^a 
depends on the two variables 



z = 



Mi 



M, 



K 



y 



'^PK-Pe 

Ml 



2Ee 

Mk' 



(2.5) 



where E^r {Eg) is the pion (positron) energy in the kaon 
rest frame, and Mk indicates the mass of the decaying 
kaon. Alternatively one may also use two of the Lorentz 

invariants 

t^{PK-P7vf, U={pK-Pef, S = (Ptt +Pe)^- (2.6) 

Then the distribution reads 

P^"Hy,z)= (y , z) I /i") (t) p , (2.7) 

^ See [3] for the general Kes form factor decomposition 
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with 



r ■ 



dydz p^°\y,z). (2.8) 



■D 



The kinematical density is given by 

A^°\y,z) =4{z + y-l)il-y)+ r,{^y + 3^-3) 



4r^ +?'e(?'7r - re), 



(2.9) 



where 



m; 



Ml 



Mi 



The boundaries of the domain of integration T) (Dahtz 
plot) in (2.8) can be found in Sect. 4. 

Virtual photon exchange as well as the contributions 
of the appropriate electromagnetic counterterms change 
the form factor [3] , 



and the distribution (2.7) has to be replaced with 
Q{y,z)=UAf\y,z)\F+{t,v)\\ 



(2.11) 



(2.12) 



The full form factor F+(i, v) depends now also on a second 
kinematical variable v as it cannot be interpreted anymore 
as the matrix element of a quark current between hadronic 
states. The variable v is taken as w = {pk — Pe)"^ for isTij 
and s = [pn + Pe)"^ for K^^. Diagrammatically, the de- 
pendence on the second variable is generated by one-loop 
graphs where a photon line connects the charged meson 
and the positron. 

The form factor Fj^{t,v) contains infrared singulari- 
ties due to low-momentum virtual photons. They can be 
regularized by introducing a small photon mass M^. The 
dependence on an infrared cutoff reflects the fact that 
F^{t,v) cannot be interpreted as an observable quantity 
but has to be combined with the contributions from real 
photon emission to arrive at an infrared- finite result. 

It is convenient to decompose Fj^{t, v) into a structure- 
dependent effective form factor f+{t) and a remaining part 
containing in particular the universal long-distance correc- 
tions [3]. To order q, the full form factor is given by 

F+{t,v) = [l + ^r{v,mlM^-M^)\ f+{t), (2.13) 

where M denotes the mass of the charged meson. Ex- 
pressed in terms of the functions Fc, A, -^2 defined in [3], 
F can be written as 

r{v,ml,M^;M^) = Fc{v,ml, M"^; M^) 



+ ri{v,ml,M^) 
+ r2{v,mlM^). 



(2.14) 



The explicit expressions for F^ Fi, F2 are displayed in 
Appendix A. The function Fc, containing a logarithmic 
dependence on the infrared regulator Mj, corresponds to 



the long-distance component of the loop amplitudes which 
generates infrared and Coulomb singularities. In the case 
of the decay, the Coulomb singularity is outside the 
physical region, while it occurs on its boundary for the K'^ 
decay. The other terms represent the remaining nonlocal 
photon loop contribution. 

Note that the effective form factor f+{t) depends only 
on the single variable t. This can be achieved [3] in the 
case of decays by the decomposition defined by (2.13) 
and (2.14). The explicit form of f^°''~{t) and f+^'^\t) 
will be reviewed in the next section. 

In order to arrive at an infrared-finite (observable) re- 
sult, also the emission of a real photon has to be taken into 
account. The radiative amplitude A4^ can be expanded in 
powers of the photon energy E^, 



(2.15) 



'(-1) ^-""(0) 
where 

(2.16) 

to the non- 
radiative amplitude M, and thus to the full form factor 

-F-f(t, u). Upon taking the square modulus and summing 
over spins, the radiative amplitude generates a correction 
p^{y,z) to the Dalitz plot density of (2.12). The observ- 
able distribution is now the sum 



Gauge invariance relates and Mj^-^ 



P{y,z) = Q{y,z) +p^{y,z). 



(2.17) 



Both terms on the right hand side of this equation contain 
infrared divergences (from virtual or real soft photons). 

Upon using (2.13) and expanding to first order in a, the 
observable density can be written in terms of a new kine- 
matical density Ai [3], and the effective form factor f+{t) 
defined in (2.13), 



p{y,z)=MSEwAi{y,z)\f+{t)f' 



(2.18) 



where we have pulled out the short-distance enhancement 
factor [5] 

Sew-.^ Sew{Mp,Mz). (2.19) 
The kinematical density Ai is given by [3] 

A,{y,z) = A^°\y,z) [l + A'^{y,z)] + A'^{y,z). (2.20) 

The function A^^{y, z) arises by combining the contribu- 
tions from lA^^'.ijP and F{v,ml, W^; M^). Although the 
individual contributions contain infrared divergences, the 
sum IS finite. The factor AY^{y,z) originates from aver- 
aging the remaining terms of [see (2.15)] and are 
infrared-finite. Note that both A^^{y, z) and AY^{y, z) are 
sensitive to the treatment of real photon emission in the 
experiment. A detailed analysis of these corrections for the 
if Jjj decay was performed in [3] . The analogous discussion 
for the ifgg case will be given in Sects. 4 and 5. 

Finally, integration over the Dalitz plot allows one to 
define the infrared-safe partial width, from which one ex- 
tracts eventually the CKM element | Vus \ ■ With the linear 
expansion of the effective form factor. 



/f "(t) = /f "(0) 1 + 



t 



(2.21) 
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the infrared-finite decay rate 



3.1 Formal expressions 



3(7) J 



r{K 7re+i/e) + r{K 7re+z^e7) (2-22) The exphcit form of /f "'^ (t) is given by [7] 



can be expressed as 



(2.23) 



where 



-fx = y dydzAi{y,z)\^l 

V 

= ao + ai A^'^ + 02 (A^'^)' 



7r± 



(2.24) 



In principle, one could easily go beyond the linear 
approximation (2.21) for the determination of the phase 
space integral. Indeed, the curvature of the form factor, 
which has been neglected in (2.21), is determined by (nu- 
merically unknown) coupling constants arising at 0{p^) 
in the chiral expansion [6]. A measurement of this curva- 
ture term in future experiments would be highly welcome. 
However, in view of the present experimental and theoret- 
ical situation, we restrict ourselves to the linear approx- 
imation (2.21). In our analysis, wo arc using the experi- 
mentally determined values of the slope parameters. This 
method [3] minimizes the uncertainties in the determina- 
tion of the phase space integrals for the time being. 

In order to extract \Vus\ at the ^1% level, we have to 



at 



provide a theoretical estimate of the form factor /_ 
t = and of the phase space integral in presence of isospin 
breaking and electromagnetic effects. We devote the next 
two sections to these tasks. 



3 The form factors f^"" (t) and /.^^'^ 



'-I- 



In this section we review the structure of the Ke3 form 
factors in the framework of chiral perturbation theory, in- 
cluding contributions of order (with isospin breaking) 
[7] and e^p^ [3], as well as effects in the isospin limit 
[6,8]. 

It is convenient [3] to write the effective form factor as 
the sum of two terms. 



(3.1) 



The first one represents the pure QCD contributions (in 
principle at any order in the chiral expansion) plus the 
electromagnetic contributions up to order e^p^ generated 
by the non-derivative Lagrangian 



(3.2) 



Diagrammatically, they arise from purely mesonic graphs. 
In the definition of ff (t), we have included also the 
electromagnetic counterterms relevant to Tr^-ry mixing. The 
second term in (3.1) represents the local effects of virtual 
photon exchange of order e^p^. 



/f (i) = 1 + ^HK+Mt) + lHK+r,{t) + Hko^- (t) 

+ V3e^^^HKn{t)-HKr,{t)] + ..., (3.3) 

where the; ellipses indicate contributions of higher orders 
in the chiral expansion (see below for the inclusion of the 
0{p^) term in the isospin limit). The function HpQ{t) 
[7,9] is reported in Appendix B. The leading order -k^-t} 
mixing angle e*-^^ is given by 



r(2) _ 



V3 



ma 



fn = {mu + rnd)/2. 



(3.4) 



The local electromagnetic term takes the form [3] 



n 



!7r2 (' 



327r2 



3 -Flog ■ 



3 log ■ 



(3.5) 



The parameter Kl^il-i) denotes the renormalized (scale de- 
pendent) part of the coupling constant K12 introduced 
in the effective Lagrangian of order e^p^ [10] describing 
the interaction of dynamical photons with hadronic de- 
grees of freedom [11,12]. The "leptonic" couphngs Ai, Xq 
have been defined in [4]. The coupling constant Ag(/x) is 
obtained from Ag(/^) after the subtraction of the short- 
distance contribution [3] , 



^6(m) 



A|D + X6-(^), 



(3.6) 



where 



p2 

e^Xf = -_ log ^ = 1 - 5ew(M„ Mz), (3.7) 

which defines [5] also the short-distance enhancement fac- 
tor SEw{Mp, Mz) to leading order. Including also leading 
QCD correction [5], it assumes the numerical value 



SwN = 1.0232. 



(3.8) 



We list here also the contributions to the form 

factor (t). Displaying only terms up to 0{p'^), the 

mesonic loop contribution is given by [3] 



(2) 



.(4) , ^(4) 



EM J 



5 1 
-HKT,{t) + -HKnit) 



+ .. 



(3.9) 



The pure QCD part of this expression was given in [7], the 
inclusion of electromagnetic contributions to the meson 
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masses and the additional contribution of 0{e^p^) due 
to TT^-T] mixing [11], were added in [12]. The sub- leading 
contributions to the TT^-r] mixing angle entering in (3.9) 
are 



-(4) 



2e(2) 



3(47rFo)2(M2 - M2) 
|(47r)2 64[3L7 + iS(M)] {M^, - M^) 



Ml 



Ml 



- M^^{Ml - Ml) log + Ml{Ml - 2,Ml) log ■ 

- 2M|,(M| - 2M2) log - 2Ml{Ml - M^) , 

A* J 

(3.10) 



and 



(4) _ 



2 V3aM| 



K 



1087r(M2-M2) 



X <! 2(4^)^ - QK;Ui) + iKl{^Ji) + 2Kl{^) + 2Kl{^i) 

(3.11) 

The local electromagnetic contribution for K'^^ is given by 



-9Z 



/:^+-°=47ra[2i^[2(/x)--Xi 
-32^(3 + l°S^^ 



3 log 



.(3.12) 



What is still missing in the expressions (3.3) and (3.9), 
is the contribution of order . Neglecting isospin break- 
ing effects at this order, the form factors of both processes 
receive an equal shift which has been calculated rather re- 
cently [8,6] in terms of loop functions (containing some of 
the Li) and certain combinations of the coupling constants 
Ci [13, 14] arising at order in the chiral expansion. For 
our purposes, we will need only the value of this contribu- 
tion at f = [6], 



/f-(0) 



= -8 


fMl-Ml\^ 


p6 





+ Aoops(m)- 



(3.13) 



3.2 Numerical estimates 

In view of the subsequent application to the extraction of 
1 14s I from partial widths, we report here numerical 
estimates for the vector form factor f^"^ at zero momen- 
tum transfer (t = 0). We recall here that in principle also 
the slope parameter A^'^ can be predicted within chiral 
perturbation theory. However, due to the relatively large 
uncertainty induced by the low energy constant Lg(Mp), 
we shall use the measured value of A^'^ in the final anal- 
ysis. 



Apart from meson masses and decay constants, which 
lead to negligible uncertainties, the vector form factor de- 
pends on a certain number of parameters (quark mass ra- 
tios and low energy constants), whose input we now sum- 
marize. 

For the quark mass ratio e^^^ defined in (3.4) we use 

[15] 

e(2) = (1.061 ± 0.083) x 10"^ (3.14) 

This number is consistent with the one obtained from a 
p^ fit [16] of the input parameters of chiral perturbation 
theory within the large errors of the latter analysis. 

For the particular combination of L, entering in (3.10), 
we take 



3L7 + Ll{Mp) = (-0.33 ± 0.08) X 10" 



(3.15) 



which is again consistent with the analysis of order p^ in 
[16]. 

For the relevant combination of electromagnetic low 
energy couplings appearing in (3.11), we use [17] 

K'^iMp) := (-6/^3 + + 2K5 + 2K6Y{Mp) 

= (5.7 ±6.3) X 10-^ (3.16) 

while for the coupling constant K12 entering in the purely 
electromagnetic part (3.5, 3.12) we take [18]: 



K{^{Mp) = (-4.0 ± 0.5) X 10" 



(3.17) 



Finally, for the (unknown) "leptonic" constants we may 
resort to the usual bounds suggested by dimensional anal- 
ysis: 



\Xl{Mp)\ < 1/(4^)2 ~ 6.3 X 10- 



(3.18) 



An alternative strategy will be discussed in Sect. 7. 

The above numerical input allows us to evaluate the 
form factor for both K'^tt~ and J^+vr" transitions. To order 
p'^, the QCD part (3.3) of the form factor at t = is 
uniquely determined in terms of physical meson masses 
(apart from a tiny contribution proportional to the leading 
order n^—rj mixing angle): 

f^"""' (0) = 0.97699 ± 0.00002. (3.19) 

Using (3.17) and (3.18), we find 

0.0046 ± 0.0001 ± 0.0008 ± 0.0003 



J + 



= 0.0046 ± 0.0008 



(3.20) 



for the local electromagnetic contribution to the form fac- 
tor. The errors given in the first line of (3.20) correspond 
to the uncertainties of KI2, Xi and Xq. In this term, the 
relative uncertainty is almost exclusively due to the poor 
present knowledge of Xi. Despite this, in the final result 
for f^°-^~ (0) this is an effect of only 0.08%. 

Combining the values given above, we obtain the result 
at 0{p'^, {md - TO„)p2, eV): 



fj^ ^ (0) = 0.9816 ± 0.0008. 



(3.21) 
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To this value, we have to add the contribution (3.13) of 
order p^, which suffers from a much larger uncertainty. Be- 
fore turning to this issue, we also list the corresponding re- 
sults for the if+TT^ form factor at 0{p^, (md— m„)p^, e^p^) 
[3]: 



f'^''''{0) = 1.0002 ±0.0022, 



= 0.0032 ± 0.0016, 
/:^+'^°(0) = 1.0034 ±0.0027. 



(3.22) 
(3.23) 
(3.24) 



3.3 The contribution 



Being the largest source of theoretical uncertainty in the 
extraction of |14s|j the p^ contribution (3.13) deserves a 
separate discussion. The loop part is given by [6] 



Aoops(Mp) = 0.0146 ± 0.0064. 



(3.25) 



The quoted error reflects the uncertainty in the p cou- 
plings L\ (contributing at order p^ through insertions in 

one- loop diagrams), as well as a conservative estimate of 
higher order effects [6] . Concerning the local contribution 
in (3.13), 



/f"(0) 



local 



Ml 



Ml 



[C[2(M,) + CJ4(M,)], 



(3.26) 

there are at present several open questions. As pointed out 
in [6] the coupUngs C[2(m) and C'^^i^^ arc experimentally 
accessible in K^-^ decays, as they are related to slope and 
curvature of the scalar form factor ]{\[t). Experimental 
efforts in this direction have started, and in the long run 
this approach will give the most reliable result. For the 
time being, following [6] we identify the estimate of short 
range contributions to /:^'^(0) given in [19] with (3.26): 



/f"(0) 



local 



= -0.016 ±0.008. 



(3.27) 



A value of this size seems to be supported by a recent cou- 
pled channels dispersive analysis of the scalar form factor 
[20], and can also be obtained by resonance saturation [21] 
for the couplings entering in (3.26), 



ares 
12 — ~ 



Fl CdC„ 



^34 - -J- 



Ml ' 



Ml 



+ *|. (3.28) 



Using [22] 

Cm = Cd = dm = F^/2V2, Mp = sf2Ms, (3.29) 

we obtain 



Inserting Ms = 1.48 GeV (scenario A of [21]), we find 



/f"(0) 



local 



= -0.012, 



(3.31) 



fully consistent with (3.27)^. 

It is important to stress here that the above methods 
do not specify the chiral renormalization scale at which the 
estimate of the relevant p^ couplings applies. This in turn 
leads to an intrinsic ambiguity in the final answer, as the 
chosen reference scale \i = Mp = 0.77 GeV is somewhat 
arbitrary. The impact of this effect can be quantified by 
studying the scale dependence of C[2+C^4 (or equivalently 
of Ziioops) with renormalization group techniques [14]. We 
find Ziioops(lGeV) = 0.0043 and Ziioops(M^) = 0.0310. We 
conclude that the present uncertainty on the p^ contribu- 
tion to /f^''(0) is at least 0.01. 

Keeping in mind the above caveats, as a net effect, 
there is a large destructive interference between the loop 
part (3.25) and the local contribution (3.27) and we arrive 
at 



/^-(O) 



-0.001 ±0.010. 



(3.32) 



Adding this number to the ones in (3.21) and (3.24), we 
obtain our final values at 0{p^, {md — m„)p^, e^p^): 



ff'' (0) = 0.981 ± 0.010, 

/:f^+'^''(o) = 1.002 ±0.010. 



(3.33) 
(3.34) 



We remark here that previous analyses [3,23] of Kes, de- 
cays and \Vus\ did not include the p^ loop contribution 
^ioops(-^p), and that further work is needed to clarify 
whether the uncertainty in (3.33) and (3.34) is a realistic 



one. 



4 Real photon radiation in K^^ 

4.1 Photon-inclusive decay distribution 

We present here in detail a possible treatment of the con- 
tribution of the real photon emission process 



KUpk) TT (p^) e+{pe) VeiPu) liP'j) 



1) I 



(4.1) 



in complete analogy with the procedure proposed in [24] 
and [3] for the analysis of the -ftT^ decay. To this end we 
define the kinematical variable [25] 



X={j)u+ P-yf = {PK -Pn-Pe 



(4.2) 



which determines the angle between the pion and positron 
momentum for given energies E,^, Eg. For the analysis of 

^ We should remark here that the estimate (3.30) is not the 
complete resonance saturation result, which actually involves 
more resonance couplings [21]. It represents, however, a well 
defined starting point and further work along these lines should 
provide the size of missing contributions and an estimate of the 
uncertainty 
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the experimental data, we suggest to accept all pion and 
positron energies within the whole K^^ Dalitz plot V given 

by 

2\/r7 <y< l + Ve-r^, 
a{y)-b{y)<z<a{y) + b{y), (4.3) 



with 



where 



a{y) = 



(2 - y) (1 + re + - y) 
2(l + re-y) 



or, equivalently, 

c{z)-d{z) <y <c{z) + d{z), 

where 

(2 -z){l + r^+re- z) 



(4.5) 



c{z) = 



2{l + r^-z) 



2(1 -+r^~ z) 

and all kinematically allowed values of the Lorentz invari- 
ant X defined in (4.2). Note that this prescription excludes 
a part of the pure -fCe37 events. The situation is best ex- 
plained by Figure 1. The dotted area refers to the K^^ 
Dalitz plot, whereas the striped region shows which part 
of the projection of the K^^^ phase space onto the (y, z) 
plane is excluded. This translates into the distribution 




Fig. 1. Dalitz plot for the three and four body final states. 



X S^'^Xpk - Pn - Pe - - P^)^\M'^f , 

pol 



(4.7) 



= M^il + r^ + Ve-y- z 



+ ■ 



yz 



+ v/(y' - 4re)(^2 _ 4^^) 



}. 



8) 



In (4.7) we have extended the integration over the whole 
range of the invariant mass of the unobserved 7 system. 
The integrals occurring in (4.7) have the general form [24] 

Im,n{P-^iP-i'-:Pi^i) '■= 

1 fd^qd^k 5^'^\P-q-k) 



2nJ q° fco (pi- k+ M^/2)"'{p2-k + M^/2y 



(4.9) 



The results for these integrals in the limit My = can be 
found in the Appendix of [24]. Using the definition (4.9), 
the radiative decay distribution (4.7) can be written as 
[25] 



p''{y,z) = 

'{y,z)Io{y,z;M^) + 



2 11/ |2| fifV- |2_^^ 



647r3 



-Cmax 



(4.10) 



where the infrared divergences are now confined to^ 
Io{y,z;Mj) = 

dx 



M2 



2pe-p-^ Il,l{Pe,P7r;PK " Ptt - Pe] M^) 
Ml Io^2{Pe,P-K;PK -Ptt - Pe;My) 
mll2,o{Pe,Pn;PK-PTr-Pe;M^) ■ (4.11) 



The explicit form of the function Jo can be found in Ap- 
pendix C. The coefficients Cm n were given in Eq. (21) of 
[25]. 

The function A^^ introduced in (2.20) can now be re- 
lated to 7o by 



A'''{y,z) = 



Io{y, z; M^) + -r{s, mi, M^; M^) 



(4.12) 

An analytic expression of the integral occurring in the 
last line of (4.10) was given in Appendix B of [26] in terms 
of the quantities V^: 



/ d^ ^ ^ (^m,n^m,n — ^ ^ ^ 



(4.13) 



^ The right-hand side of the corresponding expression (6.7) 
in [3] should be multiplied by 1/4 
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Table 1. A'i^\y,z) x 10^ for K°3 decay 



\y 


0.05 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


1.05 


6.54 


36.54 


58.54 


72.54 


78.54 


76.54 


66.54 


48.54 


22.54 


1.00 




19.54 


43.54 


59.54 


67.54 


67.54 


59.54 


43.54 


19.54 


0.95 




2.54 


28.54 


46.54 


56.54 


58.54 


52.54 


38.54 


16.54 


0.90 






13.54 


33.54 


45.54 


49.54 


45.54 


33.54 


13.54 


0.85 








20.54 


34.54 


40.54 


38.54 


28.54 


10.54 


0.80 








7.54 


23.54 


31.54 


31.54 


23.54 


7.54 


0.75 










12.54 


22.54 


24.54 


18.54 


4.54 


0.70 










1.51 


i:i.51 


17.51 


13.51 


1.54 


0.65 












4.54 


10.54 


8.54 




0.60 














3.54 


3.54 





Table 2. [Ai{y, z) - Af \y, z)] x 10^ for K% decay 



Ay 


0.05 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


1.05 


1.685 


2.071 


1.513 


0.562 


-0.523 


-1.541 


-2.295 


-2.542 


-1.864 


1.00 




2.188 


1.999 


1.243 


0.237 


-0.799 


-1.651 


-2.061 


-1.582 


0.95 




1.775 


2.024 


1.464 


0.562 


-0.432 


-1.295 


-1.757 


-1.356 


0.90 






1.844 


1.524 


0.749 


-0.180 


-1.022 


-1.501 


-1.140 


0.85 








1.476 


0.856 


0.012 


-0.792 


-1.269 


-0.925 


0.80 








1.313 


0.901 


0.162 


-0.589 


-1.049 


-0.705 


0.75 










0.883 


0.276 


-0.407 


-0.839 


-0.471 


0.70 










0.772 


0.353 


-0.243 


-0.633 


-0.200 


0.65 












0.384 


-0.097 


-0.428 




0.60 














0.031 


-0.212 





As already noticed in [3], the quantity J'j{i) given in Eq. 
(A9) of [26] (which is needed for the evaluation of V7 = U^) 
contains two mistakes: the plus-sign in the last line of (A9) 
should be replaced by a minus-sign, and 1/3™'*'^ | at the end 
of the first line of (A9) should simply read (3^^^^^. The 
function introduced in (2.20) can now be written as^ 

The expressions in (4.12) and (4.14) fully determine 
the radiativcly corrected decay density Ai{y,z) (2.20). 
In order to appreciate the effect of these universal long- 
distance corrections, we report the kinematical density 
A^i^ in the absence of clectromagnetism for several in- 
dividual points of the Dalitz plot in Table 1, while the 
corresponding radiative corrections entering in (2.20) are 
displayed in Table 2. Note that the relative size of the elec- 
tromagnetic corrections for some points (especially near 
the boundary) exceeds the average shift considerably. For 
completeness, we display a sample of numerical values for 
the kinematical densities (2.9) and (2.20) also for the K'^^ 
decay mode in Tables 3 and 4. 



* Setting ^ = in the expressions of [26] amounts to neglect 
the form factor f-{t), which is an excellent approximation in 
Ke3 modes 



4.2 Phase space integrals 

Once the function Ai{y, z) is known, the numerical coeffi- 
cients cio,i,2 entering in the phase space integral (2.24) can 
be calculated by integration over the Dalitz plot. These 
are reported in Table 5 for the K^^ mode, while the cor- 
responding results for ifij can be foimd in [3]. We recall 
once again that these numbers correspond to the specific 
prescription for the treatment of real photons described in 
the previous section: accept all pion and positron energies 
within the whole Kp^ Dalitz plot D and all kinematically 
allowed values of the Lorentz invariant x defined in (4.2). 

A full evaluation of the phase space factor Ik (2.24) re- 
quires knowledge of the slope parameter. For both modes 
we employ the measured values [1]^, 

A^"'^" = 0.0291 ± 0.0018, (4.15) 
A^'^'^" = 0.0278 ± 0.0019. (4.16) 

For i^g3 decays the final numbers 

/koL=o= 0-10372, (4.17) 
Iko = 0.10339 ± 0.00063, (4.18) 

For the A'ij mode the slope parameter given in [1] has 
received a small change compared to the PDG 2000 number 
used in [3], which amounts to a negligible difference in the 
final result 
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Table 3. A'i\y,z) x 10^ for K+s decay 







\J.\JO 


n 1 ^ 




u.oo 


n 4^ 






n 7^ 




1.05 


8.10 


38.10 


60.10 


74.10 


80.10 


78.10 


68.10 


50.10 


24.10 


1.00 




21.10 


45.10 


61.10 


69.10 


69.10 


61.10 


45.10 


21.10 


0.95 




4.10 


30.10 


48.10 


58.10 


60.10 


54.10 


40.10 


18.10 


0.90 






15.10 


35.10 


47.10 


51.10 


47.10 


35.10 


15.10 


0.85 






0.10 


22.10 


36.10 


42.10 


40.10 


30.10 


12.10 


0.80 








9.10 


25.10 


33.10 


33.10 


25.10 


9.10 


0.75 










14.10 


24.10 


26.10 


20.10 


6.10 


0.70 












15.10 


11). 10 


l.l.K) 


3.10 


0.65 












6.10 


12.10 


10.10 


0.10 


0.60 














5.10 


5.10 




0.55 
















0.10 





Table 4. [Ai{y, z) - Af> {y, z)] x 10^ for K+ decay 



Ay 


0.05 


0.15 


0.25 


0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


1.05 


1.494 


1.697 


1.174 


0.313 


-0.670 


-1.593 


-2.275 


-2.486 


-1.841 


1.00 




1.708 


1.364 


0.610 


-0.320 


-1.236 


-1.946 


-2.213 


-1.638 


0.95 




1.558 


1.378 


0.732 


-0.128 


-1.006 


-1.704 


-1.983 


-1.440 


0.90 






1.356 


0.821 


0.036 


-0.796 


-1.474 


-1.758 


-1.240 


0.85 






1.321 


0.898 


0.190 


-0.593 


-1.248 


-1.533 


-1.035 


0.80 








0.971 


0.341 


-0.392 


-1.021 


-1.305 


-0.822 


0.75 










0.490 


-0.191 


-0.794 


-1.075 


-0.597 


0.70 










0.639 


0.010 


-0.566 


-0.841 


-0.348 


0.65 












0.214 


-0.333 


-0.598 


-0.020 


0.60 














-0.094 


-0.340 




0.55 
















-0.014 





Table 5. Coefficients of the K% phase space integral 





ao 


ai 


a2 


a = 


0.09390 


0.3245 


0.4485 




0.09358 


0.3241 


0.4475 



reveal that radiative corrections effectively induce a neg- 
ative shift of 0.32% in the factor 1^0 . 
On the other hand, for K^^ one finds 

/K+|a=o = 0.10616, (4.19) 
Ik+ = 0.10482 ± 0.00067, (4.20) 

corresponding to a negative shift of 1.27% induced by the 
radiative corrections. This is essentially unchanged from 
the analysis in [3]. 



5 Extraction of \Vus\ from K^^ decays 

The CKM matrix c^lcment \Vus \ can be extracted from the 
iiTgg decay parameters by 



1/2 



1 



/r""(o) 



(5.1) 



In spite of the unsatisfactory present status of the K^^ 
data, we use them here as an illustration of the application 
of the above formula. 

• With [1] 

^(•^e3(7))fit = (7.50 ± 0.08) X lO'^ s-1 (5.2) 
and (3.33), we find 

\Vus\ = 0.2153 ± 0.0011 ± 0.0007 ± 0.0022 



= 0.2153 ±0.0026, 
where the errors correspond to 

1 Ar 



(5.3) 



\Vus 



2 r A+ /+(0) 



= |Ks|(±0.5%±0.3%±1.0%). 



(5.4) 



• A more realistic estimate of the present K^^^ uncer- 
tainty is most probably given by [1] 

^(^e3(7))average = (7-7 ± 0.5) X 10*^ S'^, (5.5) 

which implies 

IKsl = 0.2182 ± 0.0071 ± 0.0007 ± 0.0022 



0.2182 ±0.0075, 



(5.6) 
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corresponding to 

= |Ks|(±3.3%±0.3%±1.0%). (5.7) 

• Finally, combining the lifetime from the PDG with 
the preliminary photon-inclusive branching ratio from 
KLOE [27] BR(/s:4(^)) = 0.384 ± 0.002stat., we find ^ 

^(^e3(7))KLOE(prel.) = (7.43 ± 0.07) X lOS'S (5.8) 

corresponding to 

\Vus\ = 0.2143 ±0.0010 ±0.0007 ±0.0022 



0.2143 ± 0.0025. 



(5.9) 



Since the present statistical precision is comparable to 
the one of the PDG fit, we expect that the experimental 
side of the problem will improve considerably as soon as 
final results from KLOE [27] and NA48 [28] will become 
available. 



6 Extraction of \Vus\ from decays 

In this section, we update our previous analysis of the iiT^ 
decay [3] in view of the new value (3.32) for the contri- 
bution of order and the recent E865 result. All other 
parameters of the isT^ analysis in [3] remain essentially 
unchanged. Due to the inconsistency between PDG 2002 
and E865 results, we analyze them separately. 



Using the PDG-fit'^ input 



(3.93 ±0.05) X lOS" 



(6.1) 



and assuming that this number refers to the inclusive 
width of Section 4 one obtains 

\Vus\ = 0.2186 ±0.0014 ±0.0007 ±0.0023 

= 0.2186 ±0.0027. (6.2) 

• The -?^e3(7) branching ratio measured by the E865 Col- 
laboration [2], when combined with the lifetime 
from the PDG, leads to the decay width 



r{K+^^^) = (4.12 ±0.08) X W^s-\ 



(6.3) 



Note that the value BR{K^^^^^) given in [2] contains 

also events outside the K^^ Dalitz plot boundary. This 
additional 0.5% contribution has been subtracted in 
(6.3) in accordance with our prescription of the treat- 
ment of real photons. Finally, we find 

\Vus\ = 0.2238 ± 0.0022 ± 0.0007 ± 0.0023 

= 0.2238 ± 0.0033. (6.4) 

^ The systematic uncertainty in the KLOE result is not yet 

known [27] 

For K^^ the difference between "fit" and "average" is not 
sizeable 



Together with \Vud\ and jV^bj as shown in (1.1), this 
number implies 

IKdl^ + + IKbj^ - 1 = -0.0024± 0.0021, (6.5) 

in rather good agreement with a unitary mixing ma- 
trix. 

The sizeable disagreement between the result of E865 
and the PDG-fit (from old experiments) calls for further 
experimental efforts in this decay channel. 



7 Combined analysis of K^^ and X+ data 

JCij and JCgg branching fractions allow for two indepen- 
dent determinations offl^°'^ (0)-|Kts|i provided one brings 
under theoretical control isospin breaking in the ratio of 
form factors at t = 0, 



r+o :=/r"°(0)//f"'^-(0). 



(7.1) 



The standard model allows a remarkably precise predic- 
tion of this quantity. The contributions of order as well 
as the couplings Xq and K12 cancel and we are left with 
the expression 



'^-1-0 



1 



(2) 



,(4) , ^(4) ^ 
'S ^EM I 



= 1.022 ± 0.003 - 167rQ;Xi, 



(7.2) 



where the ellipses in the second line stand for isospin vio- 
lating corrections arising at 0{{md — mu)p'^ ,e^p'^) in the 
chiral expansion. We expect them to shift the result at 
most by 10^'^. Also these not yet determined contribu- 
tions have been accounted for in the error given in the 
last line of (7.2). Although no theoretical estimate of the 
coupling Xi is presently available, there is no reason why 
this low energy constant should lie outside the range sug- 
gested by naive dimensional analysis (3.18). Already such 
a rough estimate of Xi shows that r+o is confined to the 
rather narrow band 



1.017 < r% < 1.027. 



(7.3) 



We emphasize that sizeable deviations from this predicted 
range could only be understood as (i) failure of naive di- 
mensional analysis for Xx (and a dramatic one) or (ii) 
failure of chiral power counting. 

On the other hand, the ratio (7.1) is related to the 
observable 



1/2 



'+0 



(7.4) 



with the caveat that the phase space factors Ik be eval- 
uated according to the same prescription for real photons 
adopted in measuring r(K'e3(7))- Once again, it is instruc- 
tive to consider several cases: 
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Using (5.2) and (6.3), we find 



^exp _ 



= 1.062 ± 0.010 ± 0.006 ± 0.003 ± 0.003 
1.062 ±0.013, 



(7.5) 

where the errors given in the first fine refer to the ex- 
perimental uncertainties of /^(Xj^^^j), r{K^^^^-^), ^ 

and A^"'^ , respectively. The outcome is clearly in con- 
flict with the prediction (7.3) of the standard model 
and indicates indeed an inconsistency of the present 
K^^ and K^^ data. This is also illustrated by Fig. 2 
where data data from K^^ (E865) and K% (PDG-fit), 
after using r^^Q as discussed above, lead to two incon- 
sistent determinations of the product fi^°^ (0) • \ Vus\- 
Taking (5.5) instead of (5.2), the resulting numbers are 



(a) 

(b) 



r+o'' = 1-049 ± 0.010 ± 0.034 ± 0.003 ± 0.003 



= 1.049 ± 0.036. 



(7.6) 



This value is consistent with (7.3) however with a large 
error caused by the big uncertainty in (5.5). 
The inconsistency is somehow mitigated when one uses 
the present PDG-fit entries for both K^^ and K°^, lead- 
ing to 



r+o'' = 1-038 ± 0.006 ± 0.006 ± 0.003 ± 0.003 



= 1.038 ±0.010. 



(7.7) 



The present confusing status is summarized in Figure 2, 

where we plot ^ (0) -iVus \ as determined from different 
K^^ and K^^ experimental input^. The points correspond- 
ing to K'^^ have been obtained by using the central value 



for 



„th 



The overall 0.78% normalization uncertainty of 



these points is not reported in the plot. 



.225 
0.22 

0.215 
0.21 

.205 



E865 



PDG-AV 



i 

PDG-FIT 



PDG-FIT I 

KLOE*i 



Fig. 2. \Vua\ ■ f+ " (0) from Ke3 modes (see text for details). 
The KLOE result is preliminary and the quoted error is star- 
tistical only [27]. 



For the analysis of forthcoming high-precision data on 
Ke3 decays we propose the following strategy: 

* Plots of this type were first used in [29] and can be found 
also in [23,27,30] 



Check the consistency of K^^ and K^^ data by com- 
paring r'^Q with the theoretically allowed range (7.3). 
Determine the low energy constant Xi from by 
inverting (7.2), 



1.022 ± ().U03(tlioor.) - ;■; 
167rQ! 



(7.8) 



(c) 
(d) 
(e) 



(We refrain from extracting a number for Xi based on 
the present data as they are apparently inconsistent.) 
Recalculate f]^°'^ from (3.5) by using the experimen- 
tally determined parameter Xi . 

Use the new number for ^ (0) in the determination 
of |T4s| as described in Sect. 5. 

Finally, one can also use the experimentally deter- 
mined Xi to improve radiative corrections to the pion 
beta decay [31], relevant for the extraction of \Vud\ 
from this mode once the PIBETA experiment finalizes 
the analysis [32]. 



8 Conclusions 

In this work, we have studied K^^ decays using chiral per- 
turbation theory with virtual photons and leptons. This 
method allows a unified and consistent treatment of strong 
and electromagnetic contributions to the decay amplitudes 
within the standard model. We have considered strong ef- 
fects up to 0{p^) in the chiral expansion. Isospin breaking 
due to the mass diff'erence of the light quarks has been 
included up to the order {md — my)p'^ ■ Electromagnetic 
effects were taken into account up to ©(e^p^). The largest 
theoretical error is generated by the contribution of 0{p^) 
inducing a 1% uncertainty in the determination of the 
Kf.^ form factors. Additional theoretical investigation is 
needed to increase our confidence in the estimate of local 
contributions at 0{p^). 

Based on our theoretical results, we have described the 
extraction of the CKM matrix element | Vus \ from experi- 
mental decay parameters and a consistency check of ifij 
and data. 

Using the recent E865 result on the branching 
ratio, we find 



|K,s| = 0.2238 ±0.0033, 



(8.1) 



being perfectly consistent with CKM unitarity. It should 
be noted, however, that the E865 ratio differs from older 
-f^ij measurements by 2.3 a. Furthermore, the E865 result 
and the present i^'g rate as given by PDG 2002 (based 
on very old data) and by KLOE preliminary results can 
hardly be reconciled within the framework of the standard 
model. Recently-completed or ongoing experiments will 
help to clarify the situation. 

Finally a short remark on \ Vud\, the second important 
source of information for the check of CKM unitarity: the 
present number for \Vud\ is extracted from super- allowed 
Fermi transitions and neutron beta decay. In principle, 
the pionic beta decay {'Kei) provides a unique test of these 
existing determinations. This decay mode is theoretically 
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extremely clean [31] and also completely consistent with 
the present analysis of K^s decays. Using the present result 
on the TTes branching ratio from the PIBETA experiment 
[32], one finds 



\Vud\ = 0.9716 ± 0.0039, 



(8.2) 



to be compared with the current PDG value shown in 
(1.1). The final result from this experiment is expected to 
reach a precision for the pion beta decay rate of about 
0.5%. Further efforts for an improvement of the experi- 
mental accuracy of TTes would be highly desirable. 
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Appendix 



ri{v,m},M'^ 



\ogR + [A- W)T{v,m}, M^) 

,2 , 



r^iv, ml M^) = - (l - ^) [ - T{v, ml M^){1 - R) 
+ \ogR 



1 



--{Z-Y)T{v,mi,M'), (A.5) 



and 



nv,mlM^) = logX. 



B Mesonic Loop Functions 

The loop function Hpg (t) [7, 9] is given by 



HpQ{t) = 



where 



(A.6) 



(B.l) 



A Photonic Loop Functions 

The photonic loop contributions to the Kfs form factors 
depend on the charged Icpton and meson masses ml M^, 
as well as on the Mandclstam variables u = {pk —pe)'^ (for 
K^^ decays) and s = (p^- +Pe)'^ (for K^^ decays). In what 
follows we denote by v the Mandelstam variable appropri- 
ate to each decay. In order to express the loop functions 
in a compact way, it is useful to define the following inter- 
mediate variables: 



(A.1) 



2^/R 

In terms of such variables and of the dilogarithm 

dt 

U2{x) = - / '^\og{l-xt), (A.2) 
the functions contributing to F {v , ml ; M^) are given 

by [3] 

rc{v,mj,M'^;M^) = 2M'^Y C{v,ml M^) 

„, Mme / XYlogX \ , , 
+ 2 log ——^ 1 + —= . (A.3) 



Civ, ml M^) 



1 



meM 1 - X2 



-hog" X + 2 logX log(l - X^) - ^ + 1 log2 R 



LhiX^) + Li2(l - ^) + Li2(l - xVR) 
V R 



(A.4) 



18(4^(*-^^^'3) 
1 1 2SpQ — t 



12 t ApQ 

2[Ap{ii) + AQ{^i)] 



[^p(/")-^q(m)] 



(B.2) 



with 



X{x, y, z) = + + — 2{xy + xz + yz), (B.3) 



EpQ =Ml + Ml, ApQ =Ml-Ml 



Ap{ii) 



M 



(47r) 



P- log 
2 o ,,2 ' 



Mj 
1? 



and 



JpQit) 



327r2 



A 



PQ 



Ml S 



log —11 - log 
t ^^Ml Apq'^'^MI 



(B.4) 
(B.5) 



Ml 



\y\t,Ml,Ml) 



X log 



\t + \^'^{t,Ml,Ml)f-Al 



PQ 



[t-\'l^{t,Ml,Ml)f-AlQ^ 



. (B.6) 



The quantity HpQ (0) appearing in the evaluation of /+ (0) 
is given by [7] 



Ml 



^-^w=-i28;^(^^+^^)'^nM 



ho{x) = 1 + 



2x 
1 — x- 



■ log X. 



(B.7) 
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For the theoretical determination of the slope parameter 
one needs the derivative of the function HpQ{t) at t = 
given by [7] 



dHpQit) 



dt 



SFoM"^^^^^ 128.^ 



1 , MpMq 
log 




hi{x) 



x-^ - 3.T^ - 3.T + 1 



2(x-iy^ 

_^ 1 /x+ 1 



logo; 



2 V a;- 1 



(B.8) 



C The function /o(2/? z; M^) for K^^ 



The analytic result for the integral Io{y, z; My) defined in 
(4.11) is given by^ 



Io{y, z; M,) = ^ log ^ log -j^ - log 

A log log — 



+ ^ [ - Li2(??i) + Li2(l/r?i) - Li2(r72) + 02(1/7/2) 



logT(Xmax)log 



1 - r(())r(.i-,,,x) 



1 - T(a:;,nax)/r(0) 
+ Li2(r(a:max)T(0)) - Li2(T(a;max)/'r(0)) 

-Li2 (r(0)2) +7rV6 

max/ \ / 1 Pt'PiT 

+ arcosh — — arcosh — — 

4P-p(P-p^ 

+ — ^2 ' 



where 



ViPrP^y 



Pe-Pn 
P-PeP-Pn 

Pe-PTT 



(C.l) 



(C.2) 



'2ptP. P-Pe P-Pn - rrq{P-p^Y - M^{P-pef ' 

(C.3) 

P = PK-P^-Pe (C.4) 



111,2 

t{x) 



1 - 27 ± + 472 - 47 



■p^ + x/2 - sj{j>vp^ + x/2)2 - m^M2 



(C.5) 



. (C.6) 



" Note that the formula for Jq given in [25] is incorrect even 
if the Errata are taken into account 
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